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c/3 , Abstract 

C^ . 

Biological soft tissues exhibit elastic properties which can be dra- 
i-rt '. matically different from rubber- type materials (elastomers). To gain 

^ I a better understanding of the role of constitutive relationships in de- 

termining material responses under loads we compare three different 
types of instabilities (two in compression, one in extension) in hyper- 
elasticity for various forms of strain energy functions typically used 
;> ■ for elastomers and for soft tissues. Surprisingly, we find that the 

^ . strain-hardening property of soft tissues does not always stabilize the 

vQ I material. In particular we show that the stability analyses for a com- 

Tlj" ' pressed half-space and for a compressed spherical thick shell can lead 

to opposite conclusions: a soft tissue material is more stable than an 
elastomer in the former case and less stable in the latter case. 

o 

._. 1 Introduction 

j^ . Elastic materials under external loads may exhibit various responses depend- 

ing on their geometry, loads, and elastic properties. For large deformations 
or for inhomogeneous and anisotropic materials as found in biological tissues, 
these responses are best described in the theory of finite deformations [1]. In 
hyperelasticity, material properties are specified by a strain-energy function 
and there is to date a large literature on the derivation and/or fitting of 
various forms of strain-energy functions [21, f3] for rubber-type materials (re- 
ferred to as elastomers). The corresponding theory for biological soft tissues 
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is more recent and is not as well established. Nevertheless, there are a few 
standard forms of strain-energy functions used to model the elastic responses 
of soft tissues. It has long been emphasized by various authors that soft tis- 
sues with their dramatic properties under extension behave differently than 
elastomers and that in many physiological systems (such as heart, arteries, 
skin, scalp, etc.) these properties are tuned to achieve specific mechanical 
goals vital for normal function and regulation |4]. A particularly revealing 
way to understand the differences between elastomers and soft tissues is to 
push the material to its extreme by bringing it to a point where a given 
configuration becomes unstable, and to compare various instabilities for dif- 
ferent geometries. Here, we consider three prototype instabilities generated 
by external loads in an incompressible isotropic elastic body made out of 
either an elastomer or a biological soft tissue material. 

We look in turn at the instabilities generated when a half-space is com- 
pressed (Section 3.1), when a spherical membrane shell is inflated (Section 
3.2), and when a spherical shell with arbitrary thickness is compressed (Sec- 
tion 3.3). These and other types of instabilities of nonlinear elasticity have 
been reviewed in an article by Gent [5] where background literature can be 
found. For the purpose of comparison, we adopt four different strain energy 
functions which are popular in literature on elastomers and on soft tissues 
namely, the Mooney-Rivlin model, the Fung model, the Gent model, and the 
one-term Ogden model. We find that soft tissues behave differently from elas- 
tomers when it comes to stability analysis. For instance half-spaces made of 
soft tissues are stable in compression, whereas half-spaces made of elastomers 
always possess a critical stretch beyond which surface instabilities develop. 
Similar conclusions are drawn for inflation instabilities of spherical mem- 
brane shells. However, thick-walled spherical shells are found more unstable 
in compression when made of soft tissues than when made of elastomers. 
The notion of a material being more or less stable than another one used 
in this paper is in terms of the critical strains where the material becomes 
unstable and not in terms of critical stresses or external loads. The general 
conclusion is that caution must be exercised when choosing an appropriate 
model for an elastomer or for a soft tissue, because their behaviours with 
respect to instabilities are not interchangeable. The next Section recalls the 
basic underlying equations, see Ogden pLj, for instance. 



2 General set-up 

2.1 Static equilibrium 

The deformation of the material body is given by x = x(X) where X and 
X describe the material coordinates of a point in the reference configuration 
and in the current configuration, respectively. Let F = Sx/SX be the de- 
formation gradient. We consider a hyperelastic incompressible body, so that 
there exists a strain-energy function W = W{F) such that the Cauchy stress 
tensor T, specifying the stress in the body after deformation, is related to 

the deformation by 

dW 
T = F^-pl. (1) 

where p is a Lagrange multiplier associated with the internal constraint of 
incompressibility. The equation for mechanical equilibrium in the absence of 
body forces is 

div T = 0, (2) 

where div denotes the divergence operator in the current configuration. 
Equation ([2]) provides, through the constitutive relationship ([1]), a system 
of three equations for the deformation x = x(X). The boundary conditions 
are imposed by prescribing the tractions Tn at the boundary, where n is the 
outward unit vector normal to the boundary. 

2.2 Strain-energy functions 

Many different general functional forms have been proposed or derived to 
model the response of elastic materials under loads [21 El E] ■ Here, for com- 
parison purposes, we choose some typical functions that have been proposed 
to model either elastomers or soft tissues. Since we focus on the role of the 
strain energy functions in instability and not on the role of possible inhomo- 
geneities or anisotropies, we restrict our attention to homogeneous isotropic 
materials. We write the energy either in terms of the principal stretches 
Ai, A2, A3 (the square roots of the principal values of FF""") or, equivalently 
for incompressible solids, in terms of the first two principal invariants of the 
Cauchy-Green strain tensors, given by 

h = Xi + \2 + A3, I2 = A2A3 + A3A1 + Aj^Aa. (3) 



Here, we limit our investigation to a few key models that capture specific 
features and are widely used (See Table 1). The main feature of interest 
for comparison is the strain-stiffening property exhibited by many biological 
soft tissues. This can be modelled either by algebraic power dependence 
(one-term Ogden model), by exponential behaviour (as in the popular Fung 
model), or by limited chain extensibility (Gent model [3 El E])- We write 
these three models with a single parameter (z/, a, /5, respectively) such that 
the classical neo-Hookean model is obtained in the limits z/ — >^ 2, a — > 0, or 
/3 — > 0. Additionally, we also use the classical Mooney-Rivlin strain-energy 
density, often used to model elastomers; however, experimental values for the 
material parameter fi vary widely in the literature and no typical range of 
values was found. 



Name 


Definition 


soft tissues 


elastomers 


Ref. 


neo-Hookcan 


Wnh = |(/l-3) 








Mooney-Rivlin 


(/i - 3) + m{/2 - 3) 
^""" 2(1 + M) 








1-term Ogden 




v>Si 


v^Z 


[lOJIIl] 


Fung 


Wfu = ^[e"(^i-^'-l] 


3 < o < 20 




nana 


Gent 


Wgc = -^ln[l-/3(/i-3)] 


0.4 </3<3 


0.005 < /3 < 0.05 


[SEllIllIS] 



Table 1: A list of strain-energy functions. Note that the materials share the 
same infinitesimal shear modulus, which without loss of generality was taken 
equal to one. The limits /x— >0,a— >0,/3— >0,z/^2all lead to the 
neo-Hookean strain-energy. 



3 Instability 

We focus on two types of instabilities. One type is related to the notion 
of limit-point instability, which typically occurs when a balloon is infiated. 
At first, the balloon is difficult to infiate, and then it may happen that its 
radius increases dramatically and rapidly, with little or no effort to produce. 
Here balloons made of elastomers behave completely differently from balloons 
made of biological soft tissues, as Osborne [16] first observed in 1909, com- 



paring "children's toys balloons" to "hollow viscera" (dog bladders). This 
instability is investigated in Section 3.2. 

The second type of instability considered in this paper is related to the 
notion of bifurcation. Bifurcation occurs at values of the material and defor- 
mation parameters for which there exist solutions to the incremental equa- 
tions of equilibrium in the neighbourhood of a finite solution. In other words, 
the onset of instability is indicated by the existence of adjacent equilibria un- 
der the same loading. To investigate that type of instabilities, we consider 
first a finite deformation x*^°''(X) and then superimpose an incremental \V7\ 
deformation ^(^''(x) as follows 

X = X^'^+ex^'\ (4) 

where e is a small parameter. It follows that the deformation gradient is now 

F = dx/dX = (1 + eF(i)) F(°), (5) 

where F^^-* = dx^^^ /dx is expressed in the current configuration. Accordingly, 
we expand the Cauchy stress tensor in e as 

T = TW+eT(i) + 0(e2), (6) 

say, and the constitutive relationship to obtain to zeroth and first orders, 

X(o) - f(°) ^^ - «W 1 t(i) - C F(^) + f(^)f(°) ^^ - »(^) 1 (7) 

where p = p^^^ + ep^^\ C is the fourth-order tensor of instantaneous elastic 
moduli, defined by 

~ aF(o)9F(o) ' ^^ 

and the derivatives of W are evaluated on F'-^-' ; see Ogden pQ for details. Then 
the stability analysis proceeds by expanding the equation for mechanical 
equilibrium ([2]) to first order in e, that is 

div T(^) = 0. (9) 

In some cases, the geometry of the problem and the deformations considered 
are simple enough so that the condition for instability related to the existence 
of solutions for Equation iQ can be written in terms of W and the Aj. 

We now consider different geometries and our four different strain energy 
functions, to study three prototype instabilities. 



3.1 The half-space in compression 

The simplest type of bifurcation is obtained by considering an incompressible 
hyperelastic half-space with a free surface, under pure homogeneous static 
deformation with principal stretch ratios Ai, A2, A3 [T71 [18]. The correspond- 
ing instability is then assumed to correspond to the appearance of wrinkles 
on the free surface, once a critical compressive stretch ratio is reached. 

Let A2 be the stretch ratio in the direction normal to the free surface. 
A plane pre-strain is associated with deformations such that A3 = 1 (axial 
compression), whereas equi-hiaxial pre- strains are obtained for either Ai = A3 
(tangential compression) or for A2 = A3 (normal compression). It follows from 
the incompressibility condition that A1A2A3 = 1 and therefore we have 

{n = —1/2 normal compression, 
n = —1 axial compression, (10) 

n = —2 tangential compression. 

Now, the half-space is occupied with an incompressible hyperelastic ma- 
terial characterized hy W = W{Xi, A2, A3). It becomes unstable and develops 
surface instability for critical principal stretch ratios such that [191 120] 



A, 



Ai 



+ XlWii - 2AiA2W^i2 + XIW22 = 0, (11) 



where Wi = dW/dXi, Wij = d'^W/{dXidXj). 

We start with the classical elastomer modelled by the Mooney-Rivlin 
energy W^j- from Table 1. Application of the previous criterion in this case 
leads to a universal condition (independent of /i) [20j : 

^n+2 ^ 3^2n+l _ ^3n + ^3^0. (12) 

Depending on n, we obtain the classical values for the critical compression 
ratio of instability, for a half-space made of a material with the Mooney-Rivlin 
(or of course, the neo-Hookean) strain energy function. Green and Zerna [T8] 
found (Ai)cr = 0.66614 under tangential compression {n = —2); Biot found 
(Ai)cr = 0.54369 under axial compression (n = —1) and (Ai)cr = 0.44375 
under normal {n = —1/2) compression. 

Next we turn to the popular Fung strain energy for biological soft tis- 
sues. We take W = Wf^ in (ITT]) and obtain after simplification the following 
bifurcation condition: 

A^+2 _^ 3^2n+l _ ^3n ^ ^3 ^ ^a^X^ - 2A?+2" + Af +^) = 0. (13) 
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Figure 1: Critical values of the stretch ratio Ai for the instability of a Fung 
elastic half-space characterized by a stiffness-hardening parameter a. 

For n = — lorn = — 2, there is a corresponding critical value a_i = 1/2, 
or a_2 ~ 0.1644, respectively, after which the bifurcation criterion has no 
positive real root, see Fig. [H We conclude that a semi- infinite solid made 
of a Fung material, under either axial or tangential compression, is always 
stable for realistic physiological values of the parameters (say a > 1/2). For 
n = —1/2 (normal biaxial compression), the criterion has a positive real root 
for all a, which however decreases rapidly toward zero, see Fig. [U hence 
for a > 3, the half-space can be compressed by more than 97 % before the 
bifurcation criterion is met. 

Next, we consider the Gent strain energy, originally proposed for rubber 
[7j but most successfully transposed to the modelling of strain-hardening 
soft tissues. We take W = Wge in (TTT]) and obtain after simplification the 
following bifurcation condition. 



Ar + 3Af -" - \t 



I 3n— 3 



+/5(3Ar -2Ar -llAf-Af-' + gAf-' + Af-' + Af "^-3Af-^+3) = 0. 



(14) 



For n = —1, n = —2, and n = —1/2, there is a critical value /3_i ~ 0.122, 
j3^2 ~ 0.06, and /3-i/2 ~ 0.170 after which the bifurcation criterion has no 
positive real root, see Fig. [2l We conclude that a Gent elastic half-space 
under axial, tangential, or normal compression can become unstable for the 



parameters values used for elastomers (0.005 < (3 < 0.05) but is always stable 
for realistic physiological values of soft tissue parameters (0.4 < [3 <?>). 

\ 
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Figure 2: Critical values of the stretch ratio Ai for the surface instability of 
a Gent material characterized by a stiffness-hardening parameter /5. 

Finally we use the one-term Ogden model, for which the bifurcation con- 
dition reads 



Ai"+" + Ai""+' - Ai"('+") - Ai'+" + V (Ai'- 



- + Ai"'^+^) = 0. (15) 



The left hand-side of this equation is equal to —1 for Ai = and to Iv for 
Ai = 1. Therefore, it admits a real root for all positive values of v and for all 
values of n. Fig. [3] shows the dependence of the critical compressive stretch 
ratios on the material parameter v. 

3.2 The thin shell in extension 

Here we consider a spherical shell subject to an internal pressure P. Let A, 
5, R denote the inner radius, the outer radius, and the radial position of a 
material surface in the reference configuration, respectively. Let a, 6, r be the 
positions of the corresponding material points in the current configuration. In 
the case of thin shells, we look for a limit-point instability, that is conditions 
under which the curve P = P{Xa) has a maximum. When the shell is close to 
that point, a small increase in pressure can result in a large, sudden increase 
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Figure 3: Critical values of the stretch ratio Ai for the instability of an Ogden 
elastic half-space characterized by a stiffness-hardening parameter z/. 



in radius; this phenomenon is famihar to those who have blown up party 
balloons. 

Before we study the case of thin shells, we consider the radial deformation 
of a shell under uniform hydrostatic pressure (applied inside or outside). We 
use spherical coordinates, in which the radial deformation is simply r = r[R), 
with deformation gradient 

F = diag{r', r/R,r/R), (16) 

where the prime denotes differentiation with respect to R. Since the material 
is incompressible its volume is preserved and 

R^-A^ = r^- a^ R^- B^ = r^- 6^ (17) 

which leads to 

l-^' = f^(l-^') = |^(l-^')' (18) 

where A = r/R, A^ = a /A, Xf, = b/B. We denote the non- vanishing com- 
ponents of T by ti = Til {radial stress), and ^2 = ^22 = ^33 {hoop stress). 
Then the stress-strain relation ([T]) reads 

ti = XiWi -p, t2 = \2W2 - p. (19) 

The only non-vanishing equation for mechanical equilibrium ([2]) in the cur- 
rent configuration is 

dti 2, , 

^ + -(ti-t2) = 0, (20) 

and a closed equation for ti is obtained by introducing the auxiliary function 
W{X) = iy(A-2,A,A): 

^ = T^'W- (21) 

As a function of the circumferential stretch A, we have 

dti W'{X) 



dX 1-A3' 



(22) 



For a shell under internal pressure P, the boundary conditions are given by 
ti{Xa) = -P and ti(Ab) = 0. Integrating ([22]), we find [21] 

,(„.£^,,. .M P-['^,^>^- (23) 
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Recall that by (ITSIl . A;, depends on Xa, so that this latter equation is a relation 
for P as a function of A^; it can be inverted to give the displacement r = r{R) 
through (IT8D . 

Now, a limit-point instability [221 l23l 1211 [25l [261 [21] occurs when there is 
a loss of monotonicity of the function ti{Xb) as a function of Xa, that is when 
the pressure-stretch curve has a a local maximum. 

For thin shells, the analysis proceeds by considering the stress to first 
order in the small parameter 6 = {B — A) /A, measuring the thickness of the 
shell (see for instance Haughton and Ogden ^21j| or Beatty |28j). Before we 
proceed with the analysis of thin shells, it is of interest to understand the 
effect of shell thickness on the instability. To do so, we use the mean-value 
theorem and the connections (IT51) to expand ([25]) to second order in 6. We 
find 



.W'(X) S^ 
A2 ^ 2A4 



^^(^) = ^^^ + ou 



A=^-2 
A 



W'iX) - (A^ - l)W"iX) 



+ 0{6'), (24) 



where A is the position of the inner radius (see Ogden [H p. 285] for the 
first-order expansion). Since the shell wall-thickness is assumed small, this 
relation also describes the stress field at every point in the shell. 

A limit-point instability occurs for Acr such that t'^(Acr) = 0. Thus we 
first differentiate (!24|) with respect to A. Then writing t[ = at order 0(5), 
we recover the classical critical circumferential stretch for thin shells [21]: it 
is Ao (say), the smallest solution larger than one of 

^"(Ao)Ao - 2W'{Xo) = 0. (25) 

Next, to explore the dependence of the critical stretch with thickness, we 
expand Acr to first order in 6 as Acr = Aq + Xct 6 + 0(5^), say. Then writing 
t[ = at order 0(5^), and making use of (1251) . we find that Acr is given by 
the surprisingly simple equation: Aq — 2Acr Ag — 1 = 0. It follows that 

Acr = Ao + 4^^ + 0(^')- (26) 

The first order correction f[26l) shows that universally (independently of the 
constitutive relation), the critical stretch of limit-point instability increases 
with thickness for thin shells. In other words, making a spherical membrane 
slightly thicker always makes it more stable in inflation, whatever it is made 
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of. Higher-order corrections depend explicitly on the choice of W and no 
universal result is available. From now on, we focus on membrane shells and 
neglect the corrections due to S (hence, Acr is now identified with Aq given 
by(l25D). 

The limit-point instability is readily found for a neo-Hookean thin shell, 
for which W{X) = 2X^ + A~^ - 3, as [2H] Acr = 7^/^ (the neo-Hookean curve 
ti(A) is shown as the limiting case in Fig. H]). Past this critical value, the 
membrane continues stretching with reduced pressure. For certain materials, 
the pressure-stretch curve may present a maximum, followed by a minimum; 
in that case, once the maximum is reached, and the pressure is increased, 
the stretch will "jump" to a significantly higher value. This phenomenon is 
illustrated on Fig. 4 by the horizontal dotted line; it is called an inflation 
jump. Note that first, a limit-point instability is of course necessary for an 
inflation jump to occur. 

We can now investigate the existence of limit-point instability and of in- 
flation jump in strain-hardening materials. This analysis has been performed 
by various authors who noted that the limit-point instability disappears as 
the strain-hardening parameter is increased [H dUl ESI EH]- Here, we review 
and expand such results and compute the exact values of the parameters 
where such instabilities disappear. 

We start with a Mooney-Rivlin material and observe that as /i increases to 
/icr the limit-point disappears and the curve ti(A) becomes strictly increasing. 
This critical point is found by solving t[ = t'{ = 0, which gives 



/icr = ^-— ^ 0.21446, Acr = (19 + gVuV^^ - 1.84073. (27) 

5(19 + 6a/II) 1/3 , cr V ; K J 



The situation is similar for Fung materials (see Fig. Hj), where we can 
readily identify the critical values of the parameters: 

1 (92 + 12 765)^^^(3 + 765) ^ ^ 

acr = — ^ — F^ ^ 0.06685, 28) 

48 57 + 7^65 ^ ^ 



Acr = ^(92 + 12^65)^/^ ~ 1.69355. (29) 

Hence when a > 0.067, as is the case for soft biological tissues, there are no 
limit-point instabilities [2^, in accordance with Osborne's early observations 
1161. 
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(3cr = ^^ — _,/3 / ^zr - 0.05669, (30) 



For Gent materials, we find that the hmit-point instabihties disappear 
when P > Per, given by 

(10 + V93)'^'(3 + V93) 
3 315 + 33^93- (lO + ^93)^^^ (3 + ^93) 

and the corresponding circumferential stretch is 

Acr = (10 + V93)^/^ - 1.64262. (31) 

For instance, Gent [26] found limit-point instabilities (and infiation jumps) 
for infiated rubber shells with /? = 0.01 and /3 = 0.03. On the other hand, 
Horgan and Saccomandi [15] estimated that /3 ~ 0.44 for the aorta of a 21- 
year-old male and that /3 ~ 2.4 for the (stiffer) aorta of a 70-year-old male, 
and clearly, there are no limit-point instabilities in those cases (Note that the 
pressure-stretch curves for Gent materials are almost identical to the ones 
shown for the Fung energy and are not shown here.) 

The behaviour for an Ogden material is slightly different (Fig. [5]). Here 
again, the limit-point instability disappears rapidly (at Uc = 3, below any 
realistic physiological values). The asymptotic limits for ti(A) as A ^ oo are 
however different (0, 2, and cxd for < z/ < 3, z/ = 3, and u > 3 respectively). 
Note finally that there is no infiation jump for any value of u. 

We conclude that for soft biological tissues, the critical parameter values 
are far below any typical range of physiological values and the limit-point 
instability is unhkely to be observed. As noted repeatedly by Humphrey and 
co-workers [29l [301 El] , this observation should be kept in mind when a strain 
energy function is chosen in numerical simulations of soft tissues, and when 
designing artificial soft tissues for experiments. The choice of a rubber-like 
strain energy in the former case, or of an elastomer in the latter case, might 
lead to instabilities which do not actually exist in the prototype soft tissue. 

3.3 The shell under compression 

Finally we consider the case of a shell of arbitrary thickness under compres- 
sion, and analyse its stability with respect to axisymmetric perturbations in 
the usual (r, 9, ip) spherical coordinates. The stressed state x^°^ is found ex- 
plicitly from the computation of the strains and stresses in a spherical shell 
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neo-Hookean case (a = 0) 



Figure 4: Pressure- stretch curve for a Fung material. The hmit a = is 
the neo-Hookean material. When < a < «„, the system exhibits a limit- 
point instability and an inflation jump (see dotted line when a = Q;cr/2). For 
a ^ ttcr, the limit-point instability disappears. 
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Figure 5: Pressure-stretch curve for a one-term Ogden material. The neo- 
Hookean material corresponds to z/ = 2. For 2 < z/ < 3, the system exhibits 
a limit-point instability, and ti(A) ^ as A ^ oo. For a = 3, ti(A) —>■ 2 
as A — » oo. For u > 3, ti(A) -^ oo with A and the limit-point instability 
disappears. 
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done in the previous Section. Once this radial stressed state is known we 
introduce an axisymmetric perturbation x^^^ which reads 



X 



W = [u{r,e),v{r,e),Of 



(32) 



where u, v are independent of </?. The gradient F*^^^ can be exphcitly com- 
puted and the condition (jH]) further simphfied by first using the incompress- 
ibihty constraint and then expanding u, v in Legendre polynomials as 

u{r,e) = J2Un{r)Pn{cose), v{r,e) = ^ K(r)^P„(cos^), (33) 

n n 

where P„(cos6') are the Legendre polynomials, see [33] for details. After fur- 
ther simplification, a single fourth-order linear ordinary differential equation 
for f/„ can be derived 



dr \ dr^ 



c, 



dr"^ dr / 



(34) 



where 



Cs — r L\2vi 
d 
dr 



C2 = r^^/:i2i2 + 4r3£i2i2, 



{2v? + 2n- l)Ci22i + 2n(n + l)Cu22 



Co = {n + 2){n-l) 
d 



n{n + l)£iiii -{n^ + n- l)/:2222 - -^2233 - -^2121 

^ -r^>-i2i2 
dr'^ 



+ r-- (2i2i212 + -^1221 — -^2121 — -^2222 + -^2233) 

dr 

+ {71^+71 + 1)£2121 - >Cl221 - 2£i212 + -^2222 " C 



2233 



(35) 
The boundary conditions 

T^^^n = on r = a, T^^^n = -P^^^n on r = fe, (36) 
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read explicitly 

d^Un ^ (fUn r. dUn ^ ^r r. d^Un 2 dUn (n^ + n-2)^^ ^ 

ar-^ ar"^ ar otr^ r ar r^ 

(37) 

where 



d 6 

-l-'Ci212 H ' 

dr r 



D2 — ^--^1212 H -^1212, 



Di = — [-{n + n^ + A)Ci2i2 + n{n + l)Ciiii - 2n{n + l)/:ii22 + ^233 






+ {n + n^ - l)i:2222 - 2—Ci2i2r + 2n(n + l)AiW^i - \2W2 

dr 



D^ = -hn + 2){n-l) 






-^2233 + 2£i212 — -^2222 + ^"'^1212^ — A2W^2 

ar 



(38) 



The integration of Equation (!MI) takes place between r = a and r = h for 
an initial thickness A/B and the problem is to find the value of a such that 
the boundary conditions are satisfied (the outer radius 6 is a function of a). 

For this problem we use numerical techniques introduced by Haughton 
and Ogden [21] and later refined by Fu [32] and by Ben Amar and Goriely [33] , 
among others. For a Fung (exponential) strain energy, we plot the first ten 
modes at a = (neo-Hookean) , 1, 5, and 10 (strong strain-hardening effect), 
see Fig. [71 The first graph (a = 0) corresponds to Fig. [6] and has already 
been obtained and commented upon in |32i i33j where analytical expansions 
were derived for the high-number regime and for the thin-walled shell limit. 
In particular, it was shown in [33] that in the limit B /A -^ 1 the bifurcation 
curve for the mode n tends to the first positive root of 

(n + 2){n- l)/i^2 + 2{n^ + n + 7)fi^ - 3n(l + n) = 0. (39) 

The first few roots Hn are shown in Fig. [61 

Before we consider the effect of strain-hardening it is of importance to 
further comment on the neo-Hookean shell. The present graph provides 
additional information. First, it shows that the most unstable mode for thick- 
walled neo-Hookean shells is mode number 10. Second, the graph makes it 
clear that at mode 10, the critical stretch tends to a value which is higher 
than 0.66614 (the compressed half-space critical stretch value, see Section 
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3.1) as A/B -^ 0. How is this analysis compatible with the analysis of 
the neo-Hookean half-space? We expect that in the limit B —^ oo with 
A constant, the shell should be equivalent to a half-space in tangential bi- 
axial compression and that we should recover the instability discussed in 
Section 3.1. In fact, going from a thick- walled spherical shell to a half-space 
requires a double limit: not only must the shell become infinitely thick (and 
so A/B — > 0), but also the wavelength of the incremental deformation must 
be infinitesimally small compared to the thickness and radius of the sphere 
(and so n — >• oo); we checked that indeed, Acr — > 0.66614 in these limits. 

This last observation turns out to be crucial to interpret correctly the 
stabihty of a compressed Fung shell when a 7^ 0. First, the graphs in Fig. [7] 
show clearly that compressed spherical shells made of a Fung material are 
always less stable than shells made of a neo-Hookean material. This trend 
is further established analytically by determining the exact value for each 
mode in the limit A/B -^ 1 which is given by the first positive root of 

4aV^° + 2a(n2 + n + l)ii^^ + {n + 2){n- l)/i^^ - 16aV^^ 

- 2a(3n2 + 3n - 2)^^^ + 2{n^ + n + 7)fi^^ + 20aV 
+ 6a{n^ + n- l)/i^ - 3n{n + 1)^^ - 8aV - 2na{n + 1) = 0. (40) 

The analysis of these roots reveals that for thin shells the critical bifurcation 
value for each mode n increases strictly with a. In the limit of thick materials 
we see that Fung shells become unstable at stretch ratios lower than Acri = 
0.75 for a = 1, Acr5 = 0.85 for a = 5, Acrio = 0.91 for a = 10. This result 
might seem counter-intuitive in light of the analysis conducted in Section 3.1 
(see also [28]) for surface stability in compression (a neo-Hookean half-space 
is unstable when Ai < 0.66614 and a Fung half-space is always stable for 
a > 0.1644.) However, those lower bounds Acri, Acr5, Acrio correspond to 
low-mode numbers (n = 3, 2, 2 at a = 1, 5, 10, respectively), and not to the 
high-mode numbers limit necessary to reach the half-space idealization. We 
further checked that, as a increases, higher modes (say n > 15) cannot be 
excited in the limit A/B — > 0. 

We also conducted numerical investigations (not reproduced here) for the 
behaviour of shells made of other materials. The results for the Gent and 
for the one-term Ogden strain energy functions are close to those for Fung 
materials that is, a shell made of either a Gent or a one-term Ogden material 
is less stable than a thick shell made of a neo-Hookean material. Again, this 
is not contradictory with the fact that the Gent and Ogden half-spaces are 
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more stable than the neo-Hookean half-space. The analysis of the half-space 
is only relevant for high modes. 

Finally, we found that a compressed spherical shell made of a Mooney- 
Rivlin material is slightly more stable than a compressed spherical shell made 
of a neo-Hookean material (with the same asymptotic limit fl39|) when A/ B — »• 
1. Recall that the Fung, Gent, and Mooney-Rivlin materials are all stiffer in 
extension (strain-hardening effect) than the neo-Hookean material. 

4 Discussion 

One cannot help but remark that stability analysis results are difficult to 
predict in nonlinear elasticity. For instance it is by now well established that 
thin-walled spherical shells made of Fung materials are extremely stable in 
inflation. This has been proved in several contexts by Humphrey and his 
co-workers (see [30] and references therein to earlier work) to refute the hy- 
pothesis of an inflation jump instability for the development and rupture of 
intracranial aneurysms. It might also be commonly accepted that Fung ma- 
terials are extremely stable in compression, because the half-space stability 
analysis points clearly in that direction, see Section 3.1. However we demon- 
strated here that the opposite conclusion is reached for thick spherical shells 
in compression. Whereas generic strain energy functions may be suitable 
to describe some properties of materials under loads, special care should be 
taken when trying to describe instabilities in nonlinear elastic materials. 
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Figure 6: A neo-Hookean shell of inner and outer radii A and B becomes 
unstable with a critical strain A^, the largest possible values of Xa (modes 
n = 2, 3, 4, 5, 6, 8, 10, 15, 20, 30, 40, 50 are shown). The critical mode 
is the first excited mode. For instance at A/B = 0.85 (vertical arrow), the 
critical mode is ra = 4. Examples of shell deformations after the bifurcation 
are shown for n = 2, 3, 10, 15. Note that the amplitude of the mode has 
been chosen to show the structure of the solution and is not related to the 
mechanical problem at hand (the stability analysis is linear and there is no 
information on the mode amplitude or its sign.) 
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Figure 7: A shell made of a Fung material of inner and outer radii A and B 
becomes unstable with a critical strain A^, the largest possible values of \a 
(modes n = 2 to 10 are shown). 
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